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Abstract There is an analogy between the propagation of fields on a curved 
O spacetime and shallow water waves in an open channel flow. By placing a 

streamlined obstacle into an open channel flow we create a region of high ve- 
psj locity over the obstacle that can include wave horizons. Long (shallow water) 

waves propagating upstream towards this region are blocked and converted 
' Q ' into short (deep water) waves. This is the analogue of the stimulated Hawk- 

ing emission by a white hole (the time inverse of a black hole). The measure- 
^ ments of amplitudes of the converted waves demonstrate that they appear 

bJT) in pairs and are classically correlated; the spectra of the conversion process 

is described by a Boltzmann-distribution; and the Boltzmann-distribution is 
I determined by the determined by the change in flow across the white hole 

^ horizon. 
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1 Motivation 



Measurement of stimulated Hawking radiation 



There is a broad class of systems where perturbations propagate on an ef- 
fective {d + 1) dimensional spacetime geometry. In the literature this phe- 
nomenon is referred to as an analogue model. The first modern paper on ana- 
logue spacetime geometry was published in 1981 by W.G. Unruh pP, followed 
by Matt Visser ^26j in 1993. It was demonstrated that sound waves in a fluid 
flow propagate along geodesies of an acoustic spacetime metric. More gener- 
ally, for a single scalar field whose dynamics is governed by some generic 
Lagrangian C{da4>,(f>), the kinematics of small perturbations around some 
background solution, <f>{t,x) = (pQ{t,x.) + e(f)i{t,x) + ^ 4'2{t,^) + ■ ■ can be 
described by a minimally coupled free scalar field, (Z\g(0jj) — V{(j)o)) 01 = 0, 
where Z\g(0j,), a d'Alembertian operator with metric tensor 

d{daCl))d{dbCl))) ^ ' 

00 

an effective curved spacetime geometry [3 . Over the last 25 years the basic 
concept of analogue models has been transferred to many different media, 
and by now we know of a broad class of systems that possess an effective 
spacetime metric tensor as seen by linear excitations. Detailed background 
information and current developments can be found in [1]. 

Analogue models of gravity provide not only a theoretical but also an ex- 
perimental framework in which to verify predictions of classical and quantum 
field theory in curved spacetimes. For example, the first model, proposed by 
W.G. Unruh in 1981, is based on the fact that sound waves propagating on 
an inviscid and irrotational fluid flow satisfy the Klein-Gordon equation in 
an effective curved background |5T]. If the velocity of the fluid exceeds the 
velocity of sound at some closed surface, a dumb hole, i.e. an analogue of a 
black hole, forms, see Fig. [T] The presence of effective horizons opens up new 
possibilities to experimentally explore the black hole evaporation / Hawking 
radiation process. 
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Fig. 1 Acoustic iiorizons. Tlie propagation of sound waves in a convergent fluid 
flow exhibiting sub- and supersonic flow regions are depicted. The dashed (red) 
black / white lines, separating the sub- and supersonic regions, indicate the location 
of the acoustic black / white horizon. From the left to the right the flow velocity is 
speeding up and slowing down again. 
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There are several hinderances that one has to overcome before testing ana- 
fogue gravity systems in a laboratory experiment. Any experimental setup has 
to fall within the approximations made when deriving the analogy. For ex- 
ample, the main difficulty in implementing acoustic black hole (dumb hole) 
horizons, is to ensure that the waves obey the linear approximation through- 
out. Shock waves (sonic booms) occur far to readily at transitions between 
sub- and supersonic flows. In fact, we are all familiar with the sonic boom 
related to the shock wave generated by supersonic aircrafts. 

In 2002 it was argued that surface waves in an open channel flow with vary- 
ing depth are an ideal toy model for black hole experiments [TS]. Unruh's 1981 
paper raised the possibility of doing experiments with these analogues. One 
issue with Hawking's derivation is its apparent reliance on arbitrarily high 
frequencies, this phenomenon is commonly referred to as the trans-Planckian 



short wavelength cutoff, which obviates this difficulty. Thus the dependence 
of the Hawking effect on the high-frequency behavior of the theory can be 
tested in such analogue experiments [23l [22l [7] . While numerical studies indi- 
cate that the effect is independent of short-wavelength physics, experimental 
verification of this would strengthen our faith in the process. The presence of 
this effect in our physical system, which exhibits turbulence, viscosity, and 
non-linearities, would indicate the generic nature of the Hawking thermal 
process. Below we present all the necessary steps to understand and carry 
out such an experiment. 



2 Black & white hole evaporation process 

One of the most striking findings of general relativity is the prediction of 
black holes, accessible regions of no escape surrounded by an event horizon. 
In the early 70s, Hawking suggested that black holes evaporate via a quantum 
instability [TUl HO] • The study of classical and quantum fields around black 
holes shows that small classical as well as quantum field excitations are being 
amplified. In particular, a pair of field excitations at temporal frequency / 
are created, with amplitudes ay, /3f (Bogoliubov coefficients) related by. 



where gn is the surface gravity of the black hole, and a/ and f3f are positive 
and negative norm components |101 120] . Positive norm modes are emitted. 



^ The original derivation by Hawking radiation predicts that the quantum field ex- 
citations in the initial state — which are responsible for the late time radiation — 
have frequencies exponentially higher than the frequency associated with the Planck 
scale [51 [TT]. 




The dispersion relation of gravity waves creates a natural physical 




(2) 
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while negative ones are absorbed by the black hole, effectively reducing its 
mass. The surface gravity for a non-rotating black hole with a mass M is given 
by = 1-0 X 10^^/M [kg/s]. Equation ([2| is applicable for both stimulated 
and spontaneous emission, and at regimes where the quantum physics is 
dominant. A comparison of ^ with the Boltzmann-distribution allows one 
to associate a temperature T with the black hole, 

T = ^ = 1.2 . 10-12 • .gH [sK] = 6.03 • 10-«^[K] . (3) 

Here Mo is a solar mass, and the smallest observed black holes are of this 
order. Thus black hole evaporation is clearly difficult to observe directly [B]. 

The situation is not as challenging in an analogue gravity experiment, 
where one is dealing with table-top experiments that are under much better 
control and significantly easier to access. The question then arises as to how 
to collect conclusive experimental evidence to be assured one is dealing with 
analogue black hole evaporation, and not with some other classical or quan- 
tum process. As we will demonstrate below, the Hawking process exhibits 
in principle the following measurable characteristics: (i) the emission of field 
excitations is correlated; (ii) the spectra of the emission process is described 
by a Boltzmann-distribution; (iii) the Boltzmann-distribution is determined 
by the surface gravity at the effective horizon; and (iv) the emitted field ex- 
citations are stronger-than-classically / quantum correlated. In the following 
we will present an analogue gravity experiment in which we observe all clas- 
sical features of the Hawking process, i.e. (i)-(ni). We wiU later argue that it 
is not practical to look for (iv) due to the particular analogue system we are 
using. 



3 Experimental setup 

Our experiments were performed in a 6.2 m long, 0.15 m wide and 0.48 m 
deep flume (Fig. [2]), and were partly motivated by experiments in similar 
flumes [2 H?! HU . We set up a spatially varying background flow by placing 
a 1.55 m long and 0.106 m high obstacle in the flume. 

Particular care was taken to design an obstacle to minimize, or avoid, flow 
separation. Especially downstream of the obstacle as the flow slows down 
it has the tendency to separate and create a recirculating flow, see Fig. [3j 
Initially our obstacle was modeled after an airplane wing with a flat top 
and a maximum downstream slope of 5.2 degrees designed to prevent flow 
separation, with a profile given by 



H{x) = 2ai 1 — X — exp(— 6 x) 



(4) 
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where a = 0.094 of a meter and h — 5.94 per meter. However the gradual 
change in slope along the down stream side of the profile, as well as the 
absence of any sharp transitions, were not sufficient to fully prevent fiow 
separation. To address this issue we added a constant slope along the backside 
of the obstacle. The plate is 0.81 meters long. It tapers at a 4.5 degree angle 
on each end so as to create a smooth transition from obstacle to plate and 
then from the plate to the bottom of the fiume. The gradual slope eliminated 
apparent flow separation. Maximum fiow velocity occurs at the crest of the 



Fig. 3 Image of flow 
separation. The image vi- 
sualize the flow behavior 
at the lee side of an ob- 
stacle with a trapezoidal 
profile. The visualization 
technique utilizes emerged 
neutrally buoyant parti- 
cles. The motion of the 
particles during the ex- 
posure time causes streak 
lines indicating approxi- 
mately the velocity field 
of the flow. 
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obstacle. In order to reduce wave tunneling effects between the effective black 
and white hole horizons, the crest of the obstacle was extended. This was done 
by cutting the obstacle at the crest and adding a plate (15 cm in length) to 
join the sections. The extended flat section at the crest of the obstacle resulted 
in a region of relatively uniform maximum flow velocity. The final obstacle is 
displayed in Fig. |4] 





Fig. 4 Obstacle used for experiments: (la) and (lb) curved parts motivated by 
airplane wing; (2) flat aluminum plate to further reduce flow separation; and (3) flat 
top aluminum plate to reduce wave tunneling effects. 



We used particle imaging velocimetry [T] to determine the flow rate q, and 
to verify the suppression of flow separation. In this technique, small neutrally 
buoyant, tracer particles are added to the fluid. A short light pulse from each 
of two lasers (with different colors, red and green, say) is focussed into a 
narrow sheet within the fluid, the two pulses being separated in time by a 
few milliseconds. The light scattered in a direction normal to the sheet by 
the tracer particles within the sheet is focussed so that it forms an image 
of the particles in a monochrome CCD camera. Each particle produces two 
images, separated by a distance that is a measure of the component of the 
velocity in the plane of the light sheet with which the particle is moving; the 
distance between the two images is therefore a measure of the component of 
the local fluid velocity. Analysis requires that the pair of images belonging 
to a given particle be identified, and this is achieved by a cross-correlation 
technique based on the assumption that particles within a small interrogation 
area are moving with approximately the same velocity. As a result one obtains 
the flow velocity as a function of height v{h), and the flow rate is given by 

q = JI^" v{h) dh. 

Shallow water waves of approximately 2 mm amplitude were generated 2 m 
downstream of the obstacle, by a vertically oscillating mesh, which partially 
blocked the flow as it moved in and out of the water. The intake reservoir had 
flow straighteners and conditioners to dissipate turbulence, inhomogeneous 
flow, and surface waves caused by the inflow from the pump. The flume was 
transparent to allow photography through the walls, and the experimental 
area was covered to exclude exterior light. 

We are interested in the excitations propagating on the background flow in 
our setup. We measured and analyzed the variations in water surface height 
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using essentially the same techniques as in The water surface was il- 
luminated using laser-induced fluorescence, and photographed with a high- 
resolution (1080p) monochrome camera. The fluoresces served to scatter the 
light to the sides where it could be photographed, to sharply delineate the 
surface, since the mean free path of the laser in the dyed water was less than 
1 mm, and to suppressed the speckle which bedevils all laser illuminated ob- 
jects. 

The camera was set up such that the pixel size was 1.3 mm, the imaged 
area was 2 m wide and 0.3 m high, and the sampling rate was 20 Hz. The 
green (532 nm) 0.5 W laser light passed through a Powell lens to c reate a 
thin (~ 2 mm) light sheet (Fig.|5|. Rhodamine-WT dye was dissolved in the 
water, which fluoresced to create a sharp (< 0.2 mm) surface maximum in 
the light intensity. We interpolated the intensity of light between neighboring 
pixels to determine the height of the water surface to subpixel accuracy. 



Fig. 5 Surface wave 
detection. Diagram of 
light-sheet projection for 
surface wave detection: 
(1) water with dye, (2) 
Powell lens, (3) light 
sheet, and (4) fluorescing 
water surface. 




4 Quasi-particle excitations 



The excitation spectrum of gravity waves on a slowly varying background 
flow is well understood and has a dispersion relation given by, 

(l^) •tanh(2^fcM , (5) 
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with the frequency, / — 1/uj, where is the wave period; the wavenumber, 
where fc = 1/A is the wavelength; g is the gravitational acceleration, and h 
the depth of the fluid. We neglect surface tension and viscosity. We classify 
waves according to the value of 27r k h. Our waves all had wavelengths longer 
than about 2.1m (still water wavelengths), and surface tension would only 
play a role for waves with wavelengths less than about 1 cm. 

For 2tt kh < 1 the dispersion relation can be approximated by / = ^/ghk. 
These shallow water waves (called that because their wavelength 1/k is much 
longer than the depth of the water h) have both group and phase speed 
approximately equal to \/]}h. For 2'Kkh > 1, the dispersion relation is ap- 
proximated by / = {gk/2ny^^ . The group speed of these deep water waves 
is approximately half the phase speed, and both vary as the square root of 
the wavelength. For a given water depth, both the group and phase speeds of 
deep water waves are less than the group and phase speeds of shallow water 
waves. 



Fig. 6 Background noise. 
Fourier transform of water 
surface in flat bottom 
flume without waves; 
q = 0.039 m^/s and h = 
0.24 m. Fluctuations lie on 
upstream (red line) and 
downstream (green line) 
branches of the dispersion 
relation. Just visible at 
/ = ±3.1, fc = are 
the second transverse 
mode branches of the 
dispersion relation. Off 
the dispersion curves, 
the background noise 
amplitudes are less than 
0.1 mm. 
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To determine the ambient wave noise in our facility, and to check the effec- 
tiveness of our procedures, we conducted an experiment without the obstacle 
in place and with no wave generation. The space and time Fourier trans- 
form of the noise match the dispersion relation for this flow {q = 0.039 m^/s 
and h — 0.24 m) extremely well (Fig. |6|. In general, the amplitude of the 
Fourier components has a noise level of less than 0.2 mm away from the dis- 
persion curves. The apparently elevated noise energy crossing the k axis at 
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/ = ±3.1 Hz is due to the second transverse mode branch of the dispersion 
relation (the first transverse mode has a node at the location of the light 
sheet) . 

In [T^ Schtzhold and Unruh argued that the equation of motion of shal- 
low water waves can be cast into a wave equation on a curved spacetime 
background if the speed of the background flow varies. Assuming a steady, 
incompressible flow the velocity 

n{x) 

Here the two-dimensional flow rate q is fixed. The dispersion relation in the 
presence of a non-zero background velocity becomes, 

{f + vkf = (^^^ ■tanh{2TTkh) . (7) 

In Fig. [7j the dispersion relation is plotted for a flow typical of our experi- 
ments. Only the branch corresponding to waves propagating against the flow 
is plotted. For low frequencies, there are three possible waves, which we de- 
note according to wavenumber. The first, fc^, is a shallow water wave with 
both positive phase and group velocities, and corresponds to the wave that 
we generate in our experiments. The second, k'^^t^ positive phase veloc- 
ity, but negative group velocity. Both waves, fc^ and fc^j, are on the positive 
norm branch of the dispersion relation. The third, k-out, has both negative 
phase and group velocities, and it lies on the negative norm branch. In our 
experiment, generated shallow water waves move into a region where they 
are blocked by a counter-current, and converted into the other two waves. 
The goal of our experiment was the measurement of the relative amplitudes 
of the outgoing positive and negative norm modes to test the validity of ^ . 
(Further conversion from deep-water waves to capillary waves [H I16j are also 
possible but are not studied here.) The conversion from shallow water to 
deep water waves occurs where a counter-current become sufficiently strong 
to block the upstream propagation of shallow water waves [TSl [TOl [53] . It is 
this that creates the analogy with the white hole horizon in general relativity. 
That is, there is a region that the shallow water waves cannot access, just as 
light cannot enter a white hole horizon. Note that while our experiment is on 
white hole horizon analogues, it is because they are equivalent to the time 
inverse of black hole analogues that we can apply our results to the black 
hole situation. 
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Fig. 7 Conversion pro- 
cess. Dispersion relation 
for waves propagating 




against a flow typical of 
our experiments. A shal- 
low water wave, kin, sent 
upstream, is blocked by 
the flow and converted to 
a pair of deep water waves 
(ktut and k-^t) that are 
swept downstream. 
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5 Experimental procedure 

We are interested in the physics around white hole, not black hole horizon. 
In our particular analogue gravity system the two "outgoing" modes are now 
not on either side of the horizon but both come out downstream of the white 
hole horizon. In order to measure the effect of the horizon on incident waves, 
we sent shallow water waves toward the effective white hole horizon, which 
sits on the lee side of the obstacle. We conducted a series of experiments, with 
q = 0.045 m^/s and h = 0.194 m, and examined 9 different ingoing frequencies 
between 0.02 and 0.67 Hz, with corresponding still water wavelengths between 
2.1 and 69 meters, corresponding to 0.67 to 0.02 Hz frequencies. This surface 
was imaged at 20 frames per second, for about 200 s. In all cases we analyzed a 
period of time which was an exact multiple of the period of the ingoing wave, 
allowing us to carry out sharp temporal frequency filtering of the signals (i.e., 
eliminating spectral leakage). 

The analysis of the surface wave data was facilitated by introducing the 
convective derivative operator dt+v{x)dx- We redefine the spatial coordinate 
using. 



where x is the distance downstream from the right hand edge of the flat por- 
tion of the obstacle. The coordinate has dimensions of time, and its associated 
wave number has units of Hz. The convective derivative becomes dt + d^, or, 
in Fourier transform space, f + n. This is the term that enters the conserved 
norm. From equations (35), (36) and (87) of reference [IH] we find that the 
conserved norm has the form 



where A{f, k) is the t, ^ Fourier transform of the vertical displacement of 
the wave. In using this coordinate system the outgoing waves have an almost 
uniform wavelength even over the obstacle slope. 




(8) 




(9) 
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6 Data analysis and results 



We will illustrate the pair-wave creation process by presenting the results 
for fin — 0.185 Hz. In this case we analyzed images from exactly 18 cycles, 
measuring the free surface along approximately 2 m of the flow including the 
obstacle. We calculated from the wave characteristics, and after converting 
to ^-coordinates ([8| , the two-dimensional Fourier transformation as displayed 
in Fig. [s] (a) and (b). Note that the amplitudes of the Fourier transform at 
frequencies above and below ±0.185 Hz are very small, indicating that the 
noise level is small. 

As expected, there are three peaks, one corresponding to the ingoing shal- 
low water wavelength around k — 0, and the other two corresponding to 
converted deep water waves peaked near = 9.7 Hz and K~^f = —10.5 Hz. 
The former is a positive norm and the latter a negative norm outgoing wave, 
see equation 

In Fig. [9] we plot the wave characteristics (amplitude as function of t and 
^) filtered to give only the temporal 0.185 Hz band. Figures ^h) and[9]^c) are 
the characteristic plots where we further filter to include only k < — 1 Hz and 
K > 1 Hz respectively. These are the negative and positive norm outgoing 
components without the central peak of the ingoing wave (because of their 
very long wavelengths and the rapid change in wavelength as they ascend 
the slope, the incoming waves have a very broad Fourier transform). Recall, 
because we are only interested in counter-propagating waves, we defined pos- 
itive phase and group speeds as pointing to the left. As expected from the 
dispersion relationship, see Fig. [Tj the negative norm waves have negative 
phase velocity, while the positive norm waves have positive phase velocity. 
The complex structure in the characteristics of Fig. [9](a) arises because of 
the interference between the three components, the original ingoing wave, 
and the positive and negative norm outgoing waves. In Fig.[9](b), we see that 
the ingoing wave is blocked around ^ = 0, with only a small component 
penetrating into the region over the top of the obstacle ^ < 0. 



Our key results are presented in Fig. 10 Figure [10[a) shows the amplitude 
of the spatial Fourier transform at three selected ingoing frequencies. As the 
frequency increases, the ratio of the negative norm peak to positive norm 
peak decreases. Furthermore, the location of the positive norm peak moves 
slightly toward zero as the frequency increases, while the negative norm peak 
moves away from zero. This is to be expected from the location of the allowed 
spatial wavenumber from the dispersion plot, see Fig. [7] The red-dashed curve 
in Fig. 10 a) shows the Fourier transform in the adjacent temporal frequency 
bands for the sample case of 0.185 Hz. This is a representation of the noise, 
and is a factor of at least 10 lower than the signal in the 0.185 frequency 
band. 

To test whether or not the negative norm wave creation was due to non- 
linearities we repeated the runs at all frequencies with 50% larger amplitudes. 
The converted wave amplitudes did, in fact, scale linearly. 
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Fig. 8 The wave characteristics (a) shows the changes in the free surface. Notice the 
colors indicate the water level relative to the bottom of the tank, not that actual water 
heights. The double Fourier transformation (b) of the wave characteristics shows three 
excited frequency bands. The one at tj = represents the background (at k = 0) and a 
standing wave, refereed to as the undulation (to peaks at fc ~ ±10 Hz). The other two 
excited frequency bands at ±0.185 cycles/s correspond to the stimulated frequency 
bands. 



The crucial question is: Does the ratio of the negative to positive norm 
outgoing waves scale as predicted by the thermal hypothesis of equation ([2])? 
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Fig. 9 Demonstration of pair-wave conversion of an ingoing frequency of 
0.185 cycles/s: a, Fourier transform of unfiltered wave characteristic, b, Filtered wave 
characteristic containing only the ingoing frequency band, c and d, wave charac- 
teristics for filtered negative and positve norm modes. (The colours represent the 
amplitudes of the waves, see color bars.) 



This is shown to be the case in Fig. [I0[b) , where the norm ratios are plotted 
as a function of ingoing frequency. To calculate the norm of the outgoing 
waves we integrate / \ A{f, K)\'^/{f + k) dn over the peaks. In Fig. 10 ^b) the 
points represent the log of the ratios of these areas for each of the input 
frequencies we tested. The thermal hypothesis is strongly supported, with 
linear regression giving an inverse slope of 0.12 Hz and an offset close to zero. 
The slope corresponds to a temperature of T = 6 x 1Q~^^ K, and the offset 
is zero within our error bounds. 
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Fig. 10 Amplitudes and thermal spectrum, a, Absolute value of three different ingo- 
ing frequency bands, and typical noise level (red line), b, Ratio between negative and 
positive norm components in between 0.02 and 0.67cycles/s (red stars), and linear 
least-squares fit (red line). 



We see from Fig. |9[b,c) that the region of "wave blocking" where the 
ingoing wave is converted to a pair of outgoing waves, is not a phase veloc- 
ity horizon (where the phase velocity in the laboratory frame goes to zero). 
This is true even for the very lowest frequencies. The usual derivation of the 
temperature from the surface gravity relies on this conversion occurring at a 
phase velocity horizon. This makes the calculation of the surface gravity, and 
thus the predicted temperature uncertain. In our case estimates of the sur- 
face gravity give a predicted temperature of the same order as the measured 
temperature. What is important is that the conversion process does exhibit 
the thermal form predicted for the Hawking process. 
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This, together with the loss of irrotational flow near the horizon, and 
absence of a dependable theory of surface waves over an uneven bottom make 
prediction of the temperature from the fluid flow difficult. Our estimates — 
using the background flow parameters (i.e. flow rate and water height) to 
calculate gH = 1/2 9(c^ — v^)/dx — give us a value somewhere between 
about 0.08 and 0.18 Hz. What is important is that the conversion process 
does exhibit the thermal form predicted for the Hawking process. 



7 Conclusions and outlook 

We have conducted a series of experiments to verify the stimulated Hawk- 
ing process at a white hole horizon in a fluid analogue gravity system. 
These experiments demonstrate that the pair-wave creation is described by 
a Boltzmann-distribution, indicating that the thermal emission process is a 
generic phenomenon. It survives fluid-dynamical properties, such as turbu- 
lence and viscosity that, while present in our system, are not included when 
deriving the analogy. The ratio is thermal despite the different dispersion re- 
lation from that used by Hawking in his black hole derivation. This increases 
our trust in the ultraviolet independence of the effect, and our belief that 
the effect depends only on the low frequency, long wavelength aspects of the 
physics. When the thermal emission was originally discovered by Hawking, 
it was believed to be a feature peculiar to black holes. Our experiments, and 
prior numerical work [31] [7] , demonstrate that this phenomenon seems to be 
ubiquitous, and not something that relies on quantum gravity or Planck-scale 
physics. 

Black holes are linear phase-insensitive held ampliflers of a very peculiar 
kind [2S1 [12]. As mentioned in the introduction, the energy of the modes 
suffers an extremely severe de-amplification, going from frequencies and wave 
numbers far far higher than the Planck scale, to ones in the kHz regime for 
solar mass black holes. Nevertheless, when looking at the norms of the modes, 
they act just like any other amplifier. The Hawking effect is the quantum 
noise which must accompany any amplifier, but the characteristics of that 
noise are entirely determined by the amplification properties of the amplifier, 
which can of course be measured in the classical regime. This relation between 
the classical and quantum behavior was first pointed out by Einstein in his 
relation between stimulated and spontaneous emission, by Haus and Mullen 
in their characterization of quantum noise in a linear amplifier, and by many 
others [S]. In our case, the direct measurement of the quantum noise, with 
a characteristic temperature of the order of T = 6 x 10~^^ K, is of course 
impossible. A possible step forward is to study the behavior of quantum 
noise in analogue gravity systems in Bose-Einstein condensates, e.g. [l2l IT3] . 
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